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Quantum simulation of energy transport with embedded Rydberg aggregates 
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We show that an array of ultracold Rydberg atoms embedded in a laser driven background gas 
can serve as an aggregate for simulating exciton dynamics and energy transport with a controlled 
environment. Spatial disorder and decoherence introduced by the interaction with the background 
gas atoms can be controlled by the laser parameters. This allows for an almost ideal realization of a 
Haken-Reineker-Strobl type model for energy transport. Physics can be monitored using the same 
mechanism that provides control over the environment. The degree of decoherence is traced back 
to information gained on the excitation location through the monitoring, turning the setup into an 
experimentally accessible model system for studying the effects of quantum measurements on the 
dynamics of a many-body quantum system. 


PACS numbers: 82.20.Rp, 32.80.Rm, 42.50.Gy 

Introduction: Excitation transport through dipole-dipole 
interactions mm plays a prominent role in diverse phys¬ 
ical settings, including photosynthesis 0 0 , exciton 
transport through quantum-dot arrays [S], and molec¬ 
ular aggregates [BHH]. Of crucial importance is the com¬ 
petition between the fundamentally coherent transport 
mechanism and the coupling to the environment, which 
has been under intense scrutiny in the context of pho¬ 
tosynthesis (e.g. 0 IMS]) and recently experienced a 
resurge of interest (e.g. Often, clean studies of 

excitation transport are impeded by the large number of 
degrees of freedom in these systems, for example, strongly 
coupled vibrational modes liiin]- Ultracold atoms pre¬ 
pared in highly-excited Rydberg states exhibit similar 
dipolar state-changing interactions [2lH^ as found in 
organic molecules, but are considerably simpler to study. 
Due to their strong interactions and relative ease to con¬ 
trol using lasers, Rydberg atoms have been proposed as 
quantum simulators for quantum spin models imiiH] and 
electron-phonon interactions (2^]. Aggregates formed by 
networks of Rydberg atoms (Rydberg aggregates) [301 IM] 
are also ideally suited to the study of dipolar energy 
transport in an experimentally accessible system, as re¬ 
cently demonstrated |5^ . 

Here we study energy transport through a Rydberg ag¬ 
gregate embedded within an optically driven background 
gas that acts as a precisely controlled environment. This 
system extends the one recently used to observe diffu¬ 
sive excitation transport [32] by separating the aggre¬ 
gate degrees of freedom from those of the background 
gas. The background gas is electromagnetically rendered 
transparent for a probe beam. Only in the vicinity of the 
aggregate atoms, interactions disrupt this transparency, 
causing each aggregate atom to cast a shadow with ra¬ 
dius given by the interaction strength. We demonstrate 
parameters for which a larger absorption shadow is cast 
by the atom carrying an excitation, allowing us to infer 
its location. We show that the background gas simul¬ 



FIG. 1: (color online) Sketch of an embedded Rydberg aggre¬ 
gate. An assembly of several Rydberg atoms in a state | s) 
(large blue) and one in a state | p) (large orange) is linearly ar¬ 
ranged with spacing d in a background atomic gas (shades of 
green). These background atoms are then addressed with an 
EIT scheme (right panel), providing detection signals within 
radii Rc,b/p around each aggregate atom. 


taneously causes a back-action on the aggregate which 
can give rise to non-Gaussian disorder as well as site- 
dependent dephasing. The resulting excitation transfer 
dynamics can be described by a master equation similar 
to the one introduced by Haken-Reineker-Strobl (HRS) 
[53H55] to study the transition from coherent to incoher¬ 
ent transport. 

The experimental realization of a controllable HRS- 
type model will benefit the study of excitation transport 
in an open system, be it semi-conductors or light har¬ 
vesting. For the latter extensions to exciton-vibrational 
coupling and non-Markovian environments may be re¬ 
quired 0UnilSnilMlISZj • Finally, we show how decoher¬ 
ence in this system is intimately linked to the informa¬ 
tion obtained by the background gas acting as a quan¬ 
tum measurement device. In particular, despite strong 
aggregate-background interactions, decoherence vanishes 
if the background atoms do not allow one to infer the 
location of the excitation. 
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Scheme and model: The system we propose consists of a 
chain of N Rydberg atoms with spacing d forming the ag¬ 
gregate sketched in Fig. [T] Such an arrangement can be 
created by exciting Rydberg states from a trapped ultra¬ 
cold atomic gas using tightly focused laser beams [551155] . 
or by pulsed or chirped excitation in the dipole blockade 
regime gnun], which gives rise to spatially correlated 
Rydberg excitation patterns [421151] . N — 1 atoms are 
initially prepared in the state | s) = \vs) with principal 
quantum number v and angular momentum I = 0, while 
a single atom is excited to the state \p) = | r'p), with 
angular momentum 1 = 1. This \p) excitation can then 
migrate through the aggregate through resonant dipole- 
dipole exchange interactions [22l l25] . In addition, the 
aggregate is immersed in a gas of M background atoms, 
initially prepared in the electronic ground state | 5 ), the 
positions of which could be random or arranged in a 
regular fashion. These atoms are coupled by two laser 
fields from | g) via a short-lived intermediate state | e) 
(spontaneous decay rate Tp) to a third Rydberg level, 
|r) = \v's) [55H55] . Aggregate and background atoms 
could be the same or different atomic species. 

This system is governed by the many-body Lindblad 
master equation for the density matrix p {h = 1) 


p =-i\H,p]+^CiJ^p]. ( 1 ) 

Ol 


The Hamiltonian consists of three parts, H = ffagg + 
^EiT + -ffint, for the aggregate, the background gas of 
three-level atoms and van-der-Waals (vdW) interactions 
Eiiinniinni between atoms that are in a Rydberg state. 
The super-operator Lp [p] describes spontaneous decay 
of the background atom a from level | e), thus Cq [p] = 
OpO"^ — {O^Op+pO^O)/2 and the decay operator is La = 
with = [\k){k' |]a acting on atom a only 
and fc, k' G {g, e, r, s,p}. 

The aggregate atoms are labeled by Latin indices such 
as n and m. Restricted to the Hilbert space with a sin¬ 
gle excitation and setting the constant energy splitting 
between s and p to zero, we can write 



^ = 'EWnj7r„}(7r^l, ( 2 ) 

n^m n^m 


where 17 r„) = \ss..p..ss) (all aggregate atoms are in | s) 
except the n’th, which is in \p)) and Wnm = Cs/kn ~ 
Here C 3 is the dipole-dipole interaction strength 
and r„ is the position of aggregate atom n. We call 
eigenstates of ([^ excitons EH Eg. For simplicity we have 
ignored vdW interactions between aggregate atoms [65] . 

The Hamiltonian for the background gas in the rotat¬ 


ing wave approximation reads 


-^EIT = ^ 


2 ^9 ' 2 


• h.c. 


- Apdi^) - (Ap + A,)d(“) 


(3) 


where and Ap ^ are the probe and coupling Rabi 
frequencies and detunings respectively. Typically Hp <C 
Lie and Ap -|- Ac = 0 which corresponds to conditions of 
electromagnetically induced transparency (FIT) used for 
Rydberg atom detection [321 El]. 

Background atoms interact among themselves and 
with the aggregate through vdW interactions 


Hint 


a<i/3 aG{s,p},an 


For simplicity we assume isotropic interactions. To use 
the background gas as a probe for the state of the ag¬ 
gregate it is necessary that the interactions are state de¬ 
pendent. The interaction strength between two atoms 
a, n is - r„|) = Cp(a),ra/\ra - 

when they are in states |r) and |a) G {|s),|p)}. As 
concrete examples we consider | s) = | 43s), \p) = \ 43p) 
in ®’^Rb, with two choices |r) = |38s) or |r') = | 17s) 
for the upper state of the FIT ladder. The former real¬ 
izes power laws r]{a) = 6,6,4 for a = r,s,p, respectively, 
with \V<-^Pl\ > |y(’''^)|, due to a nearly resonant process: 
43p -I- 38s GG 41d -I- 38p [B5 ]j and the latter has r]{p) = 6 
and 

Excitation detection: On resonance (Ap ^ = 0), the back¬ 
ground gas becomes transparent for the probe beam de¬ 
scribed by Lip. However, close to the aggregate atoms, 
interactions = H^/( 2 Fp) destroy the trans¬ 

parency [5U ED (see also [55]). This creates an ab¬ 
sorption shadow around each aggregate atom, the radius 
Rc,a = i‘ 2 Cp(^a),raLp/Lliy/P^°''> of which depends on the 
state a G {s,p}, as sketched by blue (orange) circles in 
Fig.g Through this difference we can infer the location 
of the p-excitation. 

Effective aggregate model: To derive an effective model 
for the aggregate alone we proceed by adiabatically elim¬ 
inating the internal states of the background atoms fol¬ 
lowing the approach described in Ref. [67] ■ This is jus- 
tihed when the time scale on which background atoms 
would approach a steady state, set by the atomic de¬ 
cay rate l/Fp, is shorter than that for excitation trans¬ 
port 1/W{d) = d^/Cs (for details see [55]). The evo¬ 
lution of the reduced aggregate density matrix = 

J2nm Pnm ) ( 7 r„ |, in the case Ap^^ = 0, = 0 and 
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to leading order in flp, obeys: 
x(agg) 


* * ^„£r 


Heff = E 




E p 

ni 1 


K, 


ds*=E 




+ (Ka/^c)" 

1 


eff 


Tn)(' 


-^/r^ * + VcjVn. 


'^n ) ( "^n I j 


( 6 ) 

(7) 



where we have introduced the background-aggregate in¬ 


teraction Vna = yiil'’ +Y.m^n ■ Notc the imaginary 

contributions to . For the case Ap_c ^ 0, see [68] . 

The effective Hamiltonian describes a mean en¬ 
ergy shift of aggregate site n due to the interaction with 
the level | r) of the background atoms, weighted by the 
steady-state occupation of | r). The strength of the sec¬ 
ond term Q is set by the two-level atom photon scat¬ 
tering rate 7 eff ~ fl^/Tp within the critical radius of 
an aggregate atom. Imaginary off-diagonal terms in (§ 
arising from imaginary parts of Q can be interpreted 
as a contribution to the disorder [HS] , while real ones 
describe dephasing mechanisms. The relative contribu¬ 
tions of disorder and dephasing terms can be controlled 
by choosing Rydberg states with different interactions 
and through the FIT laser parameters. Eq. ([^ furnishes 
a Haken-Reineker-Strobl type model |3S] for excitation 
transport. All scenarios from dominant dephasing to 
dominant disorder can be realized by varying the inter¬ 
mediate state detuning Ap while keeping the two-photon 
detuning fixed: Ap + Ac « 0. In particular, for large 
Ap the contribution of dephasing can be significantly re¬ 
duced, see [SHI- 

In the following we analyze the influence of the dis¬ 
order and dephasing introduced by the background gas. 
More explicitly Eq. § reads Pnm — '^(y^kmPnk 
^^nkPkm) “t“ E^i H“ ^nrti)Pnm ^nm Pnml^i with. 


(rs) 


(na) 

eff 


= E„toid7’dr’'' 


and 


En = E„ff, 

7.™=E„(iii7’t+ \E'r'p - 2ReiLir’dr’'i)- 

then define distributions Pe(E„ - (E„)), Pe{enm - (inm)) 
and for the probability with which an individ¬ 

ual background atom a contributes to disorder and de¬ 
phasing in an ensemble average over background atom 
positions {En disorder from ®> ; ^nm disorder from 0, 
P^ dephasing). Both the width and the shape of these 
distributions can be controlled by the laser parameters 
and the background atom density. In Eig. [^ we show 
two examples: panel (a) corresponds to resonant EIT ex¬ 
citation (Ap^c = 0) and large background gas density, 
resulting in dominant dephasing and Gaussian distribu¬ 
tions, while panel (b) shows the case of finite intermediate 
state detuning and low density such that interactions be¬ 
tween probe and aggregate atoms are weaker. Detuning 
from the intermediate level reduces spontaneous decay 
and makes dephasing weaker than disorder. Remarkably, 


FIG. 2: (color online) Varying disorder and dephasing in 
quantum simulations of energy transport, (a) Histograms for 
disorder and dephasing (i) Pe, (ii) Re, (iii) P-y, using param¬ 
eters Q,p = 1.3 MHz, Dc = 30 MHz, Ap = Ac = 0, Fp = 6.1 
MHz, d — 19/im, W{d) = 0.24 MHz, p^g = 3.8 x 10^® 
m“®, G4,rp = —1032 MHz pm'*, Ce,rs = —87 MHz pm®, 
thus assuming an upper background level | r). (b) The same 
as (a), with parameters Dp = 12 MHz, Dc = 90 MHz, 
Ap = —20 MHz, Ac = 22 MHz, d = 24pm, W{d) = 0.1 MHz, 
pbg = 9.5 X 10*^ m“®, Ge.rp = —0.4 MHz pm®, Ge.rs = —0.1 
MHz pm®, thus assuming an upper background level \r'). 
(c) Effect of disorder on a single transport realization, using 
parameters as in (b). (d) Variance of the excitation location 
(x), fit by = Sp-. From bottom to top parameters as 

in (a), (b), (a) with Dp = 0. 


for low densities and weak interactions we find significant 
outliers in the atomic distance distribution that cause 
non-Gaussian disorder which can crucially modify exci¬ 
tation transport [5H]. By controlling the placement of 
individual background atoms using microstructured op¬ 
tical traps, even more exotic forms of disorder could be 
studied. 

The effects of disorder and dephasing on transport can 
be seen in Fig. (c), where we show a single realiza¬ 
tion of Eq. ([^ for A = 11 atoms immersed in a gas 
of randomly but homogeneously distributed background 
atoms. In a corresponding ensemble average, the spa¬ 
tial width of the excitation distribution over aggregate 
sites = (n^) — (n)^ carries the transport signatures 
Fig. I (d). Parametrizing we find ^ = 2 

for ballistic transport (Dp = 0), ^ = 1 for typical dif¬ 
fusive transport resulting from Fig. (a) and ^ = 0.69 
for sub-diffusive transport arising from the non-Gaussian 
disorder in Fig. 111(b). 

Imaging and measurement-induced decoherence: The de¬ 
gree of decoherence present in this system is intimately 
linked to the action of the background gas acting as a 
real-time probe of the aggregate, making it an appeal¬ 
ing model to demonstrate measurement-induced decoher¬ 
ence m- Since the background gas degrees of freedom 
have been eliminated in the effective model, we demon¬ 
strate this effect with simulations of the full master equa¬ 
tion, which also serve to verify model (§. We study 
an aggregate with A = 3, probed by two randomly dis¬ 
tributed pairs of background atoms, using a quantum- 
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FIG. 3: (color online) Link between absorption signature and 
excitation transport, (a) Difference in optical response be¬ 
tween dynamic and empty aggregate Ax{x) = Im[x — Xref] 
during transport, normalized by the two-level {\g}, |e)) re¬ 
sponse. Green lines indicate the location of aggregate atoms 
r„ and their thickness the population p„. (c) Snapshot of 

Ax(a:) at t = 3.7/rs. (b) Visualization of the correspond¬ 

ing two-dimensional signal izu, green bars indicate Pn. Here, 
Dp = 0.2 MHz, otherwise parameters as in Fig. 

jump Monte-Carlo technique [niE2]- The background 
atom pairs have a separation Ar = 0.3 /rm, yielding 
V(’'’')(Ar) = 730 GHz to include significant inter¬ 
actions between background atoms. We initially prepare 
the aggregate in state | tti ) and all background atoms in 
their ground state \g). 

Each background atom a heralds the arrival of an 
aggregate excitation through the optical susceptibility 
Xaif) = rp/HpTr(pi7eg^), the imaginary part of which 
yields the optical absorption. The average optical sus¬ 
ceptibility of the background gas x(a;) is approximated 
by spatial binning of the Xa(t) from many simulations. 
To monitor the excitation transport, one can infer the 
location of the \p) state by subtracting from x{x) a ref¬ 
erence signal Xref(a;) corresponding to the absorption of 
an inactive aggregate (chain of only | s) states) as in 
[32]. We see in Fig. that the resulting signal is di¬ 
rectly linked to the probability distribution of the exci¬ 
tation p„(t) =Tr(/5[| 7r„ )(7r„ I]). These simulations also 
show that background-background interactions are 
relatively benign for the chosen states and densities. 

The dephasing of the aggregate depends strongly on 
the position of the background atoms. In particular a 
given background atom only provides signihcant infor¬ 
mation on the excitation location if it is located in a ring 
between the two critical radii iic,s < r < Rc^p as visible 
in Fig. 1^ This is demonstrated in Fig. [^ where we place 
one background atom at a distance 6 from each site as 
shown in the top panels. For S < Rc^s background atoms 
permanently scatter a large number of photons, nonethe¬ 
less the aggregate dynamics proceeds coherently (panel 
a). In contrast, for i?c,s < S < Rc.p, despite a smaller to¬ 
tal number of scattered photons, aggregate decoherence 
is strong. 

The connection between information provided by the 



FIG. 4: (color online) Exciton transport in a continuously 
monitored embedded Rydberg aggregate. Geometries are 
shown in the top panels, (a) Site occupations Pn for n — 1,2,3 
(solid red, dashed blue, dot-dashed black) in a non-decohering 
case. Dp = 1.3 MHz, Dc = 30 MHz, d — 19pm, 5 = 0.6pm, 
other parameters as in Fig.[^(a). Golormatched crosses show 
the populations according to the effective model, Eq. (|^. 
The orange line counts the number of scattered photons 
?iphot(^) = / d^rp5^^Tr(p(t)[| e){e|]t»). (b) The same for 

a strongly decohering case with S — 1.5pm, other parameters 
as in (a). 

scattered photon and decoherence is explicit in the 
quantum-jump algorithm: The inset of Fig. (b) shows 
for a single realization how the state of the aggregate P 2 
(blue dashed) is linked to quantum jumps of the | e) pop¬ 
ulation of its probe atom (green). This link only occurs 
when the state of the background atom and the state of 
the aggregate are significantly entangled in the moment 
of spontaneous decay. Since this is not the case in panel 
(a), single trajectories (not plotted) there show no effect 
of quantum jumps on the state of the aggregate. 
Conclusions and outlook: We have shown that a Ryd¬ 
berg aggregate embedded in an optically coupled back¬ 
ground gas realizes a flexible quantum simulator of a 
Haken-Reineker-Strobl type model for energy transport. 
Site-dependent dephasing and disorder can be controlled 
through laser intensities, frequencies and background 
atomic density. Furthermore, this system could be ex¬ 
tended to study other fundamental features believed to 
be at play in photosynthetic light-harvesting: We have 
seen evidence for non-Markovian features and non-trivial 
relaxation when the time scale on which the background 
atoms reach their steady state is made comparable to 
transport time scales, a regime not discussed here. The 
analogue of internal molecular vibrations could be en¬ 
gineered as in Ref. [21]. Disorder distributions could 
be controlled even further using an additional class of 
background atoms m- All these features would extend 
the HRS type model proposed here to quantum simu¬ 
lations of light-harvesting processes in a similar spirit 
but with complementary technology to the proposals of 
Ref. EZIITIIEZ]. 

Decoherence of the aggregate arises through contin¬ 
uous monitoring of the location of the excitation, pro¬ 
viding a hands-on example of measurement-induced de- 
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coherence of a quantum state. Further applications of 
this system could be monitoring and decoherence of adi¬ 
abatic excitation transport involving external (motional) 
degrees of freedom [751 [75]. 

We gladly acknowledge fruitful discussions with Cenap 
Ates, Rick van Bijnen, Martin Garttner, Georg Gunter, 
Igor Lesanovsky, Sebastian Mobius, Thomas Pohl, Hanna 
Schempp and Matthias Weidemiiller. We acknowledge fi¬ 
nancial support by the Deutsche Forschungsgemeinschaft 
under WH141/1-1 and the EU Marie Curie Initial Train¬ 
ing Network (ITN) COHERENCE. 

SUPPLEMENTAL INFORMATION 

Steady state of an EIT system: Eor a single aggregate 
atom in state a interacting with one background atom 
[M] at a distance S, consider the Hamiltonian of that 
background atom, 

^EIT = + H.c.^ - ( 8 ) 


Adiabatic elimination of background atom excited states: 
Following |67j we now formally adiabatically eliminate 
the excited states of all background atoms to arrive at 
an evolution equation for the aggregate alone. The es¬ 
sential step is to divide the (many-body) Hilbert space 
into a space of interest and its complement. The former 
is represented by the projector 

4 = (16) 

n 

where the first part acts on the state space of the ag¬ 
gregate atoms and the second on that of the background 
atoms. We introduced |g) = \ g .. .g), which is the state 
where all background atoms are in \g). The complement 
of the space projected onto by Pg is thus formed by all 
many-body states involving any | e) or | r) state for the 
background atoms, projected onto by Pe = 1 — Pg. 

Segregating the total Hamiltonian from the main arti¬ 
cle into segments using the projection operator formalism 
m to first order in Op, we obtain: 

Hg = PgHPg = PgH.,,Pg, ( 17 ) 


The detuning A will in this case be given by the interac¬ 
tion with the aggregate atom as A = U*^’'“^((5). 

We can solve the corresponding master equation in¬ 
cluding spontaneous decay from state | e) for its steady 
state p and obtain 


ArlA^ + nt + {4A^ + nl)nl 

Pggy ) 4r2A2 + 8 A 2 D 2 + + f 22)2 ’ 

(9) 

4A2D2 

Pee{A) - ^^2^2 + 8^2^22 + + ^^2)2 ’ 

(10) 

PrrK ) 4P2A2 + + (f^2 + ^2)2 , 

(11) 

2 ADp(D 2 - 2 zrpA) 

Pgp ) 4r2A2 + 8A2D2 + (f^2 + f22)2 ’ 

(12) 

_ DpD,(D2 + - 2irpA) 

4r2A2 + 8A2D2 + (D2+D2)2’ 

(13) 

2 Anln^ 

Per{A)- 4P2A2 + 8A2D2 + (f]2 + f^2)2 ■ 

(14) 


We further define the steady state susceptibility 

X(A) = ^Im[p,e(A)]. (15) 

i Lp 

This expression can be used to describe the time- 
dependent absorption signal in Fig. 3: For W{d) Fp, 
we find that each background atom a adiabatically fol¬ 
lows the aggregate state through Xa{t) = Xadiab,a(l) = 

En X{Vna)Pn{t), where Vna = 1^™“^ is the 

overall interaction of the specific background atom a with 
the entire aggregate if the latter is in the state 17r„ ). 


H, = P,HP, 


= P. 




E(W'^- +!!■"■)-Apd(“) 


— (Ap -|- Ac)lj)^“^ -|- Hint + H., 


agg 


Pe, 


= P, 


V 

2 ^9 


Pg, 


C- = cl- 


(18) 

(19) 


After this segregation, the effective equation after adi¬ 
abatic elimination of the complement of our space of in¬ 
terest is [57] : 


dagg) ^ ^ p(-Sg)] + ^ £ [^(agg)]^ (20) 

eff 

ot 

HeS = a (77nh)^]^+ A Pg, (^1) 

= Lc.H^hC+- ( 22 ) 


Here La is the Lindblad operator introduced in the main 
article. Central to the effective equation are the non- 
Hermitian Hamiltonian Hnh = He — *Ea7lJ,LQ/2 and 
its inverse which we obtain now. 

It can be seen that if we neglect and Hagg as we 

will do from now, the Hamiltonian Hnh decomposes into 
the following block structure: 


Hnh — 




)H, 


(net) 


NH 


(23) 


where within the space spanned by | e)^ and 

I r )q, only. Similar blocks arise in C+ and La- 
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In that basis, reads explicitly: 

Tj{na) _ / —iTp/2 — Ap IIc/2 _ \ {0A\ 

"nh - -Ap - A, + j 

where Vna = + J2m^n i® Overall interac¬ 

tion of the specific background atom a with the entire 
aggregate if the latter is in the state | ). We further 

have 

= iL") = (x/r^,0). (25) 


Due to the block structure (23), we find the inverse of 


Hnh when we find the inverse of which is: 


= 


Vna -f 2 c /2 

-*rp /2 - Ap 


with Vncv — Vn.n Aj, A,- 


- \)Vna - 




(26) 


Using also (25), we obtain 


dagg) ^ ^ ^ ^ ^ 

eff 


(27) 


■^eff = XI XI I 7 r„ ) ( 7 r„ I, 

n a 

=X^eff“Xn)( 7 r„|, 

n 

nlVnai^l + WnaAp) 


rrina) _ 
^eff — 


p(nQ) _ 


D4 + 8U„„Ap02 + 4U„2jr2 + 4A2) 

2 ^I^a ^/bpUp 


2 U„a(rp - 2 iAp) - 


, (28) 


(29) 


in Eq. (21) and Eq. (22). In the limit Ap^c —> 0, we arrive 
at the expressions (Sjto (7) of the main article. 

In the main article we point out that even strongly ab¬ 
sorbing background atoms whose location however does 
not allow one to infer the excitation location do not 
contribute to aggregate decoherence. This is fully cap¬ 
tured in the model just derived, where Lindblad opera¬ 
tors L^'g in Eq. (7) of the main article for background 
atoms that are within a critical radius regardless of the 
excitation location n are proportional to a unit matrix 
and hence cause no decoherence (since O in the super¬ 
operator Cq[p] commutes with p, as O ^ lb). 

In a more explicit form of the HRS-type master equa¬ 


tion, defining IU„„ = 0, we have 

Pnm — ^ ^ ^{V^kmPnk ^^nkPkm) 


4 “ ^(T'm 4 ” ^nm)p nm 2 Pnm 5 

7.™ =Sdbrb+idrb - 2Re|tir’iir*i) 






^eff 


enm = Xliniwr A 


(no;) Y (7na)*i 


eff 


(30) 


for the matrix elements of = Y.nm P^rn\ 7 r„ )( tt^ |. 

We can interpret as diagonal disorder, 7 „m as de¬ 
phasing, and €nm as correction to the diagonal disorder, 
as explained in the next paragraph. 


Note that (20) is obtained as leading order of a pertur¬ 


bative expansion in C± and by assuming that Hg can be 
treated as “small” compared to He- Higher orders and 
corrections due to finite Hg within H^^^ can be incor¬ 
porated as described in m but have not been required 
here. 

Disorder correction Cnm- Since it depends on two aggre¬ 
gate site indices n, m in a non-trivial fashion, the inter¬ 
pretation of Cnm is not obvious at first. However, numeri¬ 
cal evaluation shows that for the parameters of Eig. 2 (a) 
in the main text, e„m can be approximated by 


<E(^. 


(na) 

eff 


-h: 


(ma) 


eff 


C* 4 , rp —0 


(31) 


The Gaussian distribution Fe in Fig. 2 (a) is reproduced 
by this expression with a standard deviation overesti¬ 
mating the correct one by ^ 10%. Thus, e„m can be 
approximated by sums of contributions arising from sin¬ 
gle aggregate sites, and consequently the term i{Em — 
En + enm)Pnm in Eq. (30) Can be cast into the form 

l{E’^- E'Jpnm, With E', = Ek- 2Ea^^iff“X4...=0- 
This is the reason we refer to Cnm as to a correction to 
the diagonal disorder. 


The approximation (31) does not hold true for all sets 
of parameters, in particular not for those of Fig. 2 (b). 
However, there is negligible. 

Dephasing in the continuum limit: The dephasing rates 


7 „m in Eq. (30) are obtained as a discrete sum over all 


background atoms a. In the limit of a continuous back¬ 
ground atom density (from here on denoted by Pbg), these 
rates can be given as a closed analytical expression, pro¬ 
vided that the interatomic distance of the aggregate sat¬ 
isfies d ^ Rc,rs, Rc,rp and that the probe and control 
fields are applied resonantly (Ap = 0 = Ac). In this 
continuum limit 7 „m = 7(1 ~ 5nm) actually becomes site 
independent. 

The dephasing 7 acquires three contributions, one de¬ 
pending solely on the rs-interaction, one depending solely 
on the rp-interaction, and one depending on both: 

7 “ 'Jrs 4- ^rp 4- 'JrSjTp- 


(32) 















7 


The first two contributions have the simple form 


while the third contribution reads 


AVrip 3 


"frs 3Y^r 

(33) 



3 

~ cos{7r/8)T/'^^^^’^P’ 

(34) 


'Yrs^rp 




{2 


{Vs + l)TpCl,^ + C4,.pC6,.. ^25rpf22^ _ 2Cl,p {2Cl,^ ^^l\C^,rs\ + C^,rsVi{TpCl^^f/^) 


+ (1-^3) (^C63,,^2010|C6..«|2 - C4,.p|C6..«P-^dTpfls 

2rp(17^^^^(17(5 ygCSC ^ ^ |f^4,rp| ^f^4,rp'^^TpI(174 y.pI + 

4" ^i,rs V ( ^ ) f ^4,rpC'6,rsf^c ~ y^^TpIc^^ 


(35) 


A similar analytical solution can be found for the on-site 
energy shifts (28), but since these also become indepen¬ 
dent of the site-index, and thus disorder vanishes, it is 
not shown here. 
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